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3 i Abstract 
• 1— I . _ 

^ , A recently [^j proposed method of continuous unitary transformations is used to 

' 

' eliminate the interaction between electrons and phonons. The differential equations 

for the couplings represent an infinitesimal formulation of a sequence of Frohlich- 
transformations |^]. The two approaches are compared. Our result will turn out to be 
less singular than Frohlich's. Furthermore the interaction between electrons belonging 
to a Cooper-pair will always be attractive in our approach. Even in the case where 
Frohlich's transformation is not defined (Frohlich actually excluded these regions from 
the transformation), we obtain an elimination of the electron-phonon interaction. This 
is due to a sufficiently slow change of the phonon energies as a function of the flow 
parameter. 



^New address: Max-Planck-Institut fiir KoUoid- und Grenzflachenforschung, Kantstrafie 55, 
14513 Teltow, Germany 



1 Introduction 



Around 40 years ago Bardeen, Cooper and Schrieffer developed their famous the- 
ory of superconductivity Essential for their success was the interpretation 
of an effective interaction between electrons of a many-particle system [Q. This 
interaction was assumed to be present in addition to a Coulomb-interaction be- 
tween the electrons. As Frohlich showed 1952 this effective electron-electron 
interaction can have its origin in the interaction between phonons and electrons. 
In this work Frohlich eliminated the electron-phonon interaction by a unitary 
transformation. By doing so he was able to describe the interaction mediated by 
the lattice as an effective electron-electron interaction. 

But Frohlich's approach contains some problems. He had to exclude certain 
regions in momentum space from the elimination since in these regions the trans- 
formation would become singular due to a vanishing energy-denominator. 

In this paper we will apply an elimination procedure recently developed by 
one of the authors 0. Instead of transforming the Hamiltonian in one step in 
this new approach the desired transformed Hamiltonian will be achieved step by 
step. Or more formally spoken instead of one unitary transformation a sequence 
of unitary transformations will be applied for diagonalization. In an infinitesimal 
formulation of this continuous transformation the renormalization of the coupling 
constants is described by the flow equations. In order to develop these differential 
equations some approximations will be necessary. By means of these transforma- 
tions new types of interactions mainly involving larger numbers of particles will 
be generated. They will be neglected after normal-ordering. Frohlich had to use 
similar approximations. 

Our approach has the following advantages: (i) The original phonon-coupling 
can be completely eliminated, even when the states connected by this interac- 
tion are degenerate. The continuous transformation is chosen in such a way that 
the transformed Hamiltonian does not contain any interactions between one elec- 
tron and the creation or annihilation of one phonon. These interactions are still 
present in Frohlich's approach due to normal ordering of generated interactions 
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not taken into account, (ii) Singularities in the induced electron-electron inter- 
action which appear in Frohlich's scheme, will either not appear in our approach 
or the divergencies will be less singular. In order to prove these properties the 
influence of the renormalization of the one-particle energies on the elimination 
of the electro n-phonon coupling has to be taken into account. By doing so it is 
possible to make statements which cannot be obtained by perturbation theory. 
Within this approach the way the couplings reach their renormalized value can 
be determined. By discussing the consequences of this asymptotic behavior it 
becomes clear that our approach is superior to Frohlich's. 

Our paper is organized as follows. In the following section a short review of 
Frohlich's approach will be given. In the Section ^ the flow equations describing 
the elimination of the electron-phonon interaction will be derived. In the next 
section the above mentioned approximations made by Frohlich will be applied to 
the system of differential equations. The renormalized values will be calculated. 
Also the differences in the electron-electron interactions will be discussed. Section 
1^ contains the analysis of the behavior of our transformation in the asymptotic 
regime. The fundamental differential-equation for the phonon-energies will be 
discussed. In the next section the consequences of the asymptotic behavior of 
the couplings will be given. It will be shown that the electron-phonon coupling 
always is eliminated even in the case of degeneracies. The last section contains a 
short summary and outlook. 

2 Frohlich's Transformation 

In the following sections we will often refer to the above mentioned work of 
Frohlich [^. Thus for convenience a short summary of his results will be given 



here. 



The Hamiltonian of the model will be written as 



H 




(1) 
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Here and in the following k stands for k = {k, cr}, i.e. the spin is conserved by the 
electron-phonon interaction thus no spin-subscript is needed, a*^^-* are bosonic cre- 
ation respectively annihilation operators, c''^^ denote the corresponding fermionic 
operators. : ... : denotes normal-ordering and ii^ is a constant energy. Further Mq 
is the coupling between electrons and phonons. Following the approach of Bloch 
[0 or Nordheim it is independent of the electron momentum. If there is need 
to specify Ek or Uq a quadratic dispersion for electrons and a linear dispersion 
for phonons will be assumed. Finally it should be emphasized that neither the 
Coulomb repulsion nor umklapp-processes will be taken into consideration. 

Frohlich eliminated the electron-phonon interaction in Eq. (|l|) up to order 
|Mg(0)p by expanding a unitary transformation with the Baker-Hausdorff for- 
mula 

= e-^ife^ = H+[H,S] + ^[[H, S], S] + ... 

Frohlich's ansatz 

S--=-T.M, (- ^— + ; \ -s) (2) 

assures that the relation 

i/e-p + [i^o, S\ = (3) 
is fulfilled. If degeneracies can be excluded the transformed Hamiltonian becomes 

Q k 5 

+ Kk',s '■ cl+sci'-sCk'Ck ■■ 

k,k',5 

+ irrelevant Terms , 

where the irrelevant terms are either of order \Mq\^ or represent other interactions 
not taken into account in Eq. (|I]). The occurring coefficients are given by: 

e'k = e,-Y\Mq\'(nq- + K + 1)- ^-——) (4) 

q \ Ek+q — tk — UJq Ek+q — E^ + Uq J 

< = ujq-Y.m-'-nJ + \ \ (5) 

y. X^k+q ~ UJq Ek+q — Ek — i^q / 
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yL', = m'. — — (6) 

= E + Y.nk{e^ -Sk) -Y.nknk+qVk^k+q,q (7) 

k k,q 

In the last equations the convention has been introduced that q denotes a phonon- 
wavevector, thus Uq is a bosonic occupation number whereas Uk and rik+q denote 
the fermionic ones. 



3 The Flow Equations 

In order to eliminate the electron-phonon coupling by flow equations one divides 
H 

H = H'^ + H\ 
into the phonon-number conserving part 

q k 5 

+ Yl ^k,k',S ■ cl^s4'-5Ck'Ck ■■ +E 
k,k',5 

= HP'' + H"" + H^-^ + E 
and the phonon-number violating part 

k,q 

All occurring coefficients have to be regarded as functions of the flow parameter 
I. The initial conditions of the introduced new quantities are 

M,,q{l = 0) = M*^^^_^{1 = 0) = Mq{0) ^ Mq 

Vk,k'Ai-o) = 0- 

The generator of the continuous unitary transformation will be chosen as 

V ■= Yi^k,qak,qalq + Ml+q,-qPk,qaq)cl+qCk, (8) 
k,q 
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where 

Oik,q = ^k+q ~ + h^q = ^k+q — ^k ~ ^q- 

By calculating the commutator [r], H] one obtains for dH/dl besides terms of 

the type already presented in H"^ contributions of the type W"^^^, W^~'^^^ and 
p^2e-2ph_ 'pj^gy a^j^g Qf ^]^g form 

~ E : ^V'/ ^ ^ : (9) 

^2e-2r.h ^ ^ (at_, + a,) : cI^.^+^,cl,_^,Cfc,c, : (10) 

k,k' ,q,q' 

w^'"' ~ EK^^ + ^-'^s)- (11) 

The interactions W'^"'^^^ and W'^'^~'^^^ are not included in the original model 
Hamiltonian. They describe two-phonon-processes and the interaction of a phonon 
with two electrons. Although a more realistic Hamiltonian should contain such 
interactions, we will neglect them as Frohlich did. 

The interaction (|Tl|) describes the generation and annihilation of two phonons. 
It can be transformed away by introducing 

fj = 7] + ri^'^\ 

= V + J2^<i (4^-9 ~ "-i^^-i 

9 

where fig and are real and fig{l = 0) = 0. 

These additional terms modify the equations for Ug, E, M^^q. Besides this one 
obtains the additional flow equation 

^ = -'^^qiq {Mk,gMk + g-gak + g-g " M^-g^gM k-gtt k-g) ■ 



By choosing 



= J— E {Mk,qMk+g-gak + g-g " Mk-g^gMk-gOk, 

^^9 k 



it is assured that the interaction ([IT]) will not be generated, i.e. one obtains for 
all flow-parameters / G Rq 

fiq{l) = 0. 



Therefore the desired transformed Hamiltonian becomes 

H'^ioo) = ^a;g(oo) : aja^ : +^(£^(00) - 2^nk+5Vk,k+s,s{^)) ■ ■ 
q k S 

+ Vk,k',5ioo) : 4+5cJ.,_5Cfe/Cfc : +£^(00) 

k,k',S 

+ irrelevant terms . (12) 
The renormahzation of the coefficients is described by the flow equations: 

dMu a 9 , , 



dl 



,k+q+S,sMk+S,qOik+S,q ' {iT-k+q+S — IT-k+s) 

S 

— '2Mk,qOik,q ■ X! (^*;+5^fe,fe+5,<5 ~ ^fe+g+5^fe+9,fe+9+(5,5) 

<5 

+2 ■ E Vk,k'+q,qMk>,qak',q ■ {Uk'+q — Uk' ) 
k' 

M* 

k±q^ ^ M,., gMk'+q-q/3k',q ■ {Uk'+q " Uk') (13) 



^1 k' 



dVk,k',q 
dl 
djjJn 



dl 

dsk 
'df 



Mk,qM;,_^Jk',-q - Ml^^^_^Mk>,-qak>^-q (14) 

2 • X) \Mk,q\^ak,q ■ {Uk+q - Uk) (15) 
k 

- Y,{'2nq\Mk+q,-q\^l3k,q + 2(n, + l)|Mfc,/afc,,) (16) 



dE ^ dSk v-^ Vk,k+q,q 

The generated new interactions have been neglected. Thus the last equation does 
not describe correctly the renormahzation of the ground state energy. Hence it 
will not be of further interest. 



4 Comparison with Frohlich's Method I: The 
Non-Degenerate Case 

In this first approach it will be assumed that afc,q 7^ holds for all /c, g. This could 
be realized by a finite system with appropriately chosen electron and phonon- 
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dispersion. In this case it is assured that the Frohhch-transformation is wellde- 
fined. 

In this section the flow equations will be compared with Frohlich's approach 
of Section ^. His results were only exact up to order |Mg(0)p. Thus in the flow 
equation approach terms of order |Mq(0)|^ and higher might also be neglected. 



Because of the initial condition Vk^k',q{^) = Eq. (H) shows that Vk^k',q{^) is of 



the order |Mq(0)p. Thus the lines two to four in Eq. ([131) become irrelevant. 
The same holds for the last line of this equation, i.e. is also of the order |Mg(0)|^. 
Within the order given by this approximation Eq. (|T^ can be solved exactly 

Mk,q{l) = Mg(0)e-("*+^(°)-"*(°)+"'(°))'-' + C(|Mg(0)H. (18) 

Using this solution the equations for Vk^k',q{l), ^kiP) and ujq{l) are easily inte- 
grated. Thus their renormalized values can be obtained. While efc(oo) = and 
ujq{oo) = ojq holds the interaction becomes 

T4,.'.(oo) = \MM? (-T%^ - ) ■ ^''^ 

\^k,q^ Pk'-q Pk,q ~^ '^k' ,-q / 

Actually the last equation does not represent Frohlich's result. To illustrate this 
a look at the interaction between the electrons of a Cooper-pair will be taken 
{k' = -k): 

In this case Frohlich proposes (actually independently of this specialization) 

yk,k',q = ykU,q = \Mqm'T—^^^T^—-i , (20) 



whereas the flow equations yield 



^.,-.,,(00) = -WM'i ^- (21) 

Thus a remarkable difference between continuous unitary transformations and 
the Frohlich transformation has arisen. Before discussing the origin of this dif- 
ference it should be mentioned that the difference between the two interactions 
is independent of Mg(0), i.e. holds even for arbitrarily weak coupling. 
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But it is easily seen that both approaches yield the same result for real pro- 
cesses. In this case the one-particle energy is conserved, i.e. 

£k+q — £k = £k' — £k'~q- (22) 

For these processes the interaction becomes 

V,',,^^ = \M,iO)\'- iff e,+, -e, = sy - sy^,. (23) 

K^k+q — fcfcj — 

By introducing 

B := {(A;, k' , q) : Sk+q - Ek = Sk' - Sk'-q}, 
the whole interaction can be written as 

H'-' = ( E yk,k',q + E yik'J ■ cUqci.qCk^Ck : . (24) 

\(k,k',q)£B (k,k',q)<^B J 

Both Frohlich's and the flow equation approach yield the same first term, i.e. 
both Eqs. (|^) and (|20|) yield the result ( p3D for Vj^f., q. Only the second term 
Vj[f,,^g in ( p^D differs. Thus the differences between the interactions (0) and (pO|) 
can be seen as a different generalization from processes with {k, k', q) E B to 
processes with (fc, k', q) ^ B. The accordance in V^'^fc/ ^ grants the independence 
of the self-energy part of the one-particle energies. 

Now the origin of the difference between the two interactions (|TS]) and (^) 
has to be investigated. For that it is useful to formulate the transformation 
described by the flow equation as an /-independent transformation. To do so the 
flow equations 




= [v. H] 



will be interpreted as an infinitesimal formulation of the unitary transformation 

H{1) = U{l)H{0)U\l), 

where 

U{1) := Tze/o^'C')^''. 
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Here /-ordering is defined in the same way as time-ordering. 
Thus to get rid of the /-dependence one can introduce 



e"^ := f/(oo). 
Then for 5* the following expansion holds 

/•oo 1 roo rl 

S = Si + S2 + ... = - dlr]{l) - - / dl dl'[r]{l),7]{l')] + ... (25) 

JO 2, Jo Jo 

and the neglected terms are again of order |Mg(0)|^. Here t] denotes again the 
choice (P). Using the solution (|18]) the first term of the last equation becomes the 
generator of the Frohlich-transformation, while for the second term one obtains 

k,k',q 

J <^k,q 1 f3k,q _^ 1 I 



./3fc',-g(«L + Oik,qPk',-q Oik',-q{PL + al, ) Pk,qOik',- 



^'^k+q^^k'^k'-q^k'-q- 

+ terms of the structure o)ac^c. (26) 



To derive the result (|T^) by using this modified Frohlich-transformation up to 
order |Mg(0)p only [I]fc^fcc[cfc, S-2\ has to be taken into account. Again the gen- 
erated interaction of the form (^ will be neglected^. It is easily verified that 

E^fc4cfc, '^2] = E (^fc,fc',g(oo) - V5t',g)4+g4/_5Cfc/Cfc (27) 

k k,k',q 

holds. 

Thus the change in the effective electron-electron interaction is caused by a 
change of Frohlich's generator. This modification is produced by carrying out the 
/-ordering. Thus it is a consequence of the /-dependence of rj respectively of the 
/-dependence of the coefficients of H. Or more rigorously spoken it is the flow of 
the couplings which changes the unitary transformation. 



^Actually these terms together with ~ dlrj'^'^^l) guarantee that H{oo) does not contain 



any interactions of the form (IT 
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By comparing (|I^ and it is obvious that the flow equations yield a result 
which is less singular, i.e. Vk^k',q diverges^] only for the special case ak,q = Pk'-q = 
0. Nevertheless it is not clear which transformation is better with respect to the 
neglected interactions. Although the modification 5*2 given above eliminates most 
of the singularities of Frohlich's result it is not clear by now if this modification 
5*2 is the best one can choose. By changing the coefficient of S2 independent of 
any r] a wide variety of possible corrections to V^^^./ ^ can be achieved. It is still 
an open question what the side-effect of these modifications will be. But one 
should be aware that modifications of 5*2 are equivalent to modifications of the 
the unitary transformation. Only for real processes (where energy is conserved), 
one does not have a choice to modify the effective interaction in this perturbative 
approach. 



5 The Asymptotic Behavior 

By now the elimination of the electron-phonon coupling has been considered up 
to second order. A main advantage of the flow-equations is that in certain cases 
one can go beyond perturbation theory. Then the continuous change of the single 
particle energies has the interesting effect, that for processes which are real for 
a certain value of the flow parameter /, this energy-conservation does no longer 
hold for other values of I. As Kehrein, Mielke and Neu have shown for the 
spin-boson-problem this allows a complete elimination of the coupling even if 
energy-conservation holds asymptotically, i.e. for I approaching infinity. The 
renormalization of the electron-energies Sk will be neglected, assuming that the 
electron-phonon coupling is not strong enough to cause a significant change of 
Ek- This holds at least for weak electron-phonon coupling. But special cases 
■^We will return to this point in Section ^. 

''in deriving Eqs. ( p^ and ( |l7| ) this assumption already has been used. There the l- 
dependence of the fermi-functions Uk has be neglected. This can be done for small differences 
Vfe,fc',g(oo) — Vfc^fc'.g(O). This assumption is consistent with the BCS-theory, where one neglects 
the self-energy. 
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like heavy fermions are excluded. Furthermore it is assumed that the influence 
of the lines two to flve in Eq. (0) is negligible^. Thus the dominating part of 
Mk^q becomes 

Mfc,,(/) = M,(0)e-/o'^''K^-^^+"«(''))' (28) 

It should be emphasized that this ansatz corresponds to a highly non-perturbative 
behavior. By taking into account the change of Ug in the order of |Mq(0) p the Eq. 
( |28D and hence the Eqs. (|T^) — (|l^ contain terms of any order of |Mg(0)p. This 
corresponds in the sense of perturbation theory to a summation over a particular 
subclass of diagrams. 

With these approximations only the behavior of Ug for large / has to be ex- 
amined. The key to the asymptotic regime lies in the interpretation of Eq. ([TED . 
The coupling Mk^g decays exponentially as long as ak,g does not lie in the vicinity 
of a resonance g ~ 0. Thus the behavior of Ug for large / is determined only 
by the contributions coming from small ak^g- 

To develop a more formal understanding of this point one replaces in Eq. (p!5|) 
sums by integrals and introduces the quantity T := Air ■ j^^- 

Specializing to three spatial dimensions and the case of quadratic electron- 
dispersion Ek = Y~ S^^S 



dujq 

Hi 



-r ■ \Mg 



„ „ —aq o rl 



da{kl - a') ■ + eq + ujq)e-'Io dV{-oqim^s,^.,f 

kp m 

+ da{kl - a') -i^ + Sg- ^ )e-2/o 

J-kp m J 



(29) 



where kp stands for the fermi-momentum. By introducing ujg := j Jq dl'ujg{l') one 
can use 



exp ^—2 J dl'{—aq/'m + eq±ujq)'^^ 



ujg{r)dl' ] ■exp(-2{-aq/m + eg±uJg) -l 



^It can be seen by a self-consistent analysis that this assumption does not modify the results 
of this section, i.e. the asymptotic behavior of ujq. 
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This allows to express the integral over a in terms of elementary functions and 
the error integral. For large / and for values of q which are not too large the 
integrals can be extended from — oo to +00. By neglecting the exponentially 
decaying terms the result is 

= -2 ■ r ■ |Mg(0)|2exp (-2 ujl{l')dl' + j (j^ujg{l')dl'^ 



TYl^ I 77 1 



dl 



It is worth mentioning that in leading order and besides a constant Eq. (|30|) is 
identical with the integro-differential equation describing the asymptotic behavior 
of the flow equations of the spin-boson problem 0. Following the approach of 
Kehrein, Mielke and Neu one can assume an asymptotic expansion of cUg, i.e. uiq 
decays algebraically 

UJq{l)^UJq{00)+Cq-l~'', Z/>0, / > 1, (31) 

where Cq is independent of /. Using uJq — ujq ~ the integrals in Eq. (^) 
are easily performed. The only consistent solution is the case with z/ = 1/2 and 
Cq = ±1/2. Numerical studies suggest uj^ <0, thus the ansatz (|3l|) yields 

Uqil) ^ Uqioo) + / > 1. (32) 

This solution implies that the approach to the g-dependent limit is independent 
of q. We note however that u!q{l) can be rescaled u!q{l) = q ■ Cgilq^) and Mg(0) = 
ygM(O). Then Eq. (0) yields 

dcAz) 



dz 



-2T ■ |M(0)|2 ■ exp(-2 / cl{z)dz + 2z ■ cl) 

Jo 



xm\l — ■ [cs{cs - Cs) + ^] 



with the asymptotic behavior 

1 



Cs[z) ^ cJoo + 



2^/^^ 

Assuming that this asymptotic behavior holds for z ^ zi then Eq. (|32|) holds for 
I ^ li = Zi/q^. Thus for small q it is only reached for very large /. 
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To get a better understanding of the behaviour in the asymptotic regime a 
new approach has to be chosen. Therefore one has to go back to Eq. (|30|). If uOq 
decays weaker than than only the term ~ ^^(oo) {pq—uOq) has to be considered. 
Hence the other terms will be neglected. An a posteriori justification of this 
approximation will be given later. The remaining integro-differential equation can 
be transformed into a nonlinear differential equation. One simply differentiates 
the first term of (^) and expresses the new generated functions in terms of 
l,u!q,u!'q. Introducing the function 

V{1) := / UJq{l')dl' — I ■ UJq{l), 
<j 

the result of this procedure becomes 

/2 (v" ■ V - v'"^) + h-v-v' + 2v^ -v' = (33) 

where v' = Thus the integro-differential equation (|30| ) has been replaced by 
a much more handable nonlinear differential-equation. 

To make a first connection to the ansatz (^) one sees that (|3^) corresponds 
to V = Thus the algebraic decay given above is a singular^ solution of the 
differential equation (P3|). 

In order to solve the differential equation (^^ one can make use of the fact 
that this equation is equidimensional. Thus the ansatz 

vil) = Vl-v{0. e = log///o (34) 

yields the equation 

// / /2 I o 3 / I 4 2 n i'or^\ 

V ' V ~ ~ V + 2,1] T] + T] ~ = U, (35) 

which only implicitly depends on /. Here is 77' = As a next step one can 
substitute 

77' = —777(0;), X = — - log |?7|, (36) 

which yields for 7 7^ 

^singular has to be understood as being unable to fulfill arbitrary initial conditions. 
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with 7' = Interpreting this as a differential equation for the inverse function 
one gets for 7 7^ 1/2 

dx 27 



(38) 



d-f 27-1 
By integrating the last equation one obtains 

x = 7 + ^log|27- 1| +C, C:= const. (39) 

Thus the applied transformations yield a first integral of the differential equation 
(|33D . Although the remaining Eq. (|3^) can only be solved in special cases it is 
possible to get an analytical understanding of the behavior of Uq in the asymptotic 
regime. One only has to use the fact that 7 obeys the autonomous equation (P^) 
and that the first integral (|39|) is known. The further discussion is organized as 
follows. First the solution 7 will be classified by phase space arguments. The 
properties of 7 determine the behavior of 77. Thus the behavior of the general 
solution will become clear and a connection with the singular solution will be 
made. Finally the remaining Eq. (|36D will be solved for some special cases. 

If the function j : D —>■ R is a. solution of a first order autonomous equation 
then it has exactly one of the following properties: 

(i) 7 is constant and D = R. 

(ii) 7 is injectiv and regular, i.e. 7' 7^ on D. 

(iii) 7 is regular periodic, i.e. D = R and 7' 7^ on D. 
For a proof of this Theorem see for example B. 



While the singular solution ( p2D belongs to type (i) the solution for arbitrary 
initial conditions has to be classified now. For that the function x will be inter- 
preted as function of 7, as shown in Fig. |l]. It follows from the definition (^) 
that X > Xmin '■= x{\ri\ = 1/2) (see also Fig. H). Thus the value 77 = 1/2 corre- 
sponds to the minimum of the transformation x = x{ri). Because there exists a 
lower bound Xmin the image of the function 7 = 7(x) is also bounded. Therefore 

7 cannot be bijectivQ. 

^In the following wc will sometimes refer to the inverse function of 7, which has to be 
understood as local inverse function. 



14 



0.7 - 




-0.8 -0.6 -0.4 -0.2 O 02 04 

Figure 1: x as function of 7 to different initial conditions. 

For general initial conditions the function x = x{'~f) will be periodic, as can be 
easily seen by phase space arguments. Because of x = r]"^ — log 1 7^1/2 the same 
holds for T] = ?7(log///o)- 

Knowing that in general the function 77 will be periodic some more properties 
of the solution can be derived by using Fig. |T|. First it should be mentioned that 
the singularities of the differential-equation (^) split the domain of x(7) in two 
sections. The assumption u'^ < implies 7 < 1/2, as can be easily see from Eq. 
(|34D . Thus 7 = 1/2 is the upper bound of the domain of x(7). 

Furthermore the point (7 = 0,x = Xmin) in the phase space corresponds to 
a stable trajectory. This follows from the Eqs. (^) and (0), which imply 
r^N = 0, Vn. Therefore 77 = const., which means that also x and 7 are con- 
stant. This stable trajectory is the singular solution (|32D . 

If general initial conditions x(7 = 0) > Xmin have to be fulfilled the relation 
x'(7 = 0) = still holds. But the higher derivatives of x no longer vanish. 
Therefore a small deviation from the initial condition ^(7 = 0) 
the solution 7 = 7(x) to change from the behavior of type (i) to type (iii). The 
trajectory in the phase space is rather a curve than a point. This curve connects 
the two turning-points 7;, 7.^ G {7 : x(7) = Xmin}- It corresponds to an oscillation 
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Figure 2: x as function of rj. 



of X around Xmin, i-e. of rj around 1/2. For the turning-points the relation 



holds. This implies that an oscillating trajectory cannot turn into a constant one. 
To characterize the oscillation the following consideration is useful. In Fig. |l] a 
trajectory x to the initial conditions x(7 = 0) = Ci can be seen. This phase 
space curve corresponds to an oscillation of rj along the curve x(?7) in Fig. ^ 
Therefore the turning- points 7/, 7,. in the phase space correspond to the minimum 
at ?7 = 1/2. While the point x(7 = 0) = Ci represents the greatest elongation. 
Defining ri± := x~^{Ci), where + denotes the right branch and — left left branch 
of x{ri) the following relation holds: 



Thus the amplitudes for the elongations rj > 1/2 are greater than those for 
rj < 1/2. On the other hand the 'time' for which the solution stays in the regime 
with 7] < 1/2 is longer than that for rj > 1/2. This becomes clear by \ri'\ = rj ■ I7I 
and the fact that the image of 7 is independent of the branch. But \ri\ is smaller in 
the left branch than in the right one. Therefore |?7'| is smaller in the left branch. 



dx 




1 1 



^+-2>2-^- 
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Figure 3: rj as function of log Z//o to the initial conditions 
r/(ei) = 0.6, 77' (^i) = 0. 

Altogether the deformation of the oscillations in the Figs. can be understood 
by now. 



The above mentioned differences in amplitude and 'remain-time' in the cases 
1] > 1/2 and rj < 1/2 suggest that the relation (r/^ — 1/4) =77^ — 1/4 = holds. 
Here (...) denotes an average over one period. In analogy to the prove of the virial 
theorem one can try to find a periodic function /(//, r]') fulfilling the relation 

1 



df 2 



(40) 



The condition (|40D becomes easier to handle if again a change to the variables x 
and 7 is made. Letting /(x, 7) = f{ri,ri'), one has to fulfill the condition 

dx 27 

It is sufficient to take only analytical / into account: 



<97 I 27^ 



(41) 



fix, 7) = /o(x) + ■ 7 + -Ux) ■ 7^ + ... 



Therefore one obtains the recursional relations 

/i = 

nfn - -Jn+l + ^—/j^ = 



for n > 1. 



(42) 

(43) 
(44) 
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Figure 5: 77 as function of log///o to the initial condi- 
tions r){^i) — 1.8, rj'{^i) — 0. The sequence of figures 
corresponds to a proceeding deviation from the singular 
solution. 
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The function fo{x) can be chosen arbitrarily. The ansatz /o(a;) ~ x"^, m > 2 
imphes that fn is a polynom of degree (m — 1) in x for n > 2. In order to keep 
the analysis as simple as possible only the two following cases will be discussed: 
l.case: /n(a:) = C = const.: It follows easily by induction that fn+i{x) = 2" ■ n! 
holds for n > 0. Together with Eq. (H2|) this yields 



/(x,7) = C^-^log(l-27). (45) 

Because of 7 < 1/2 the last equation is well defined. 

2.case: fn(x) = —x: Here fnix) = holds for n > 2, and therefore 

/(a;,7) = 7 - X. (46) 

Because of Eq. ( pOj) the Eqs. (^5]) and (|46|) are actually identical. Thus the 
oscillating behavior of the function rj is completely determined by the first integral 
(|39D. Simultaneously it is shown that (r^^) = 1/4. 
From the definition of v and Eq. (R3) one obtains 



From this one can conclude 



UJg{l) = CJg(oo) + Cq^, I > 1. (47) 



V/' 



Here c = Cq obeys the differential equation 



with the solution 



c(0 = exp(^/2) |"(| + V) ■ exp(-e/2)de 



Since rj is periodic in ^, this also holds for c. One notes, that rj/ 2 + 7]' = 77(1/2 — 7) 
is positive, since both factors on the right hand side have this property. Thus ujg 
decays monotonically as and the term (4/)~^ in Eq. (^) is asymptotically 
neghgible. Thus the approximations leading to Eq. ( PBD are indeed justified. 
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For the following it will be of importance to show that the average of (? with 
respect to ^ equals 1/4. From the differential equation (^) one obtains 



Since the right hand side is a total derivative the averages of (? and r/^ coincide. 



Nearly harmonic case: Here it is assumed that i] deviates only weakly from 
the singular solution (corresponding to Fig. By expanding the deviation of 
r/^ from 1/4 in powers of the amplitude A with A ^ \ the Eq. (|35|) yields the 
solution 



Thus stronger deviations from the singular solution correspond to smaller fre- 
quencies. Observe, that the term proportional to does not contain a constant. 
The same applies for the term proportional to A^ which is not given here. This 
was a first indication for us, that the average of 77^ might be 1/4 independent of 
the amplitude A. 

Instead of averaging the solution the algebraic decay of lo^ can be characterized 
in first order in A by the decay-exponents v = ^, ^ ± i^. Hence harmonic 
oscillations correspond to complex decay-exponents. 

Strongly non-harmonic case: Here the solution is characterized by a strong 
deviation from the singular solution (corresponding to Fig. |^). For the case 
^ 1/2 Eq. (^) can be expanded in powers of e if 77 is assumed to be of order 
e. In this case the solution is 



-T]^ = {c + r]){c -r]) = 2(c' + t]'){c + 7]) 



d{c + r]f 



Thus (c2) = 1/4. 

Finally two special solutions of Eq. ( |53| ) will be given. 




y 



^ + ^.{l-A^ + 0{A')} 
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where C must be positive to assure that rj has a positive slope. 
The case rj ^ 1/2 is not so easy to handle. Here t] can be expanded at its 
maximum M = r/^^^^ — |log?7max ~ Vmax- ^^e vicinity of this maximum 
7 = y/M'^ — rf holds as can be seen by expanding Eq. (^9|) . Thus a good 
approximation for this regime is 

2Me-^(«-^o) 



1 + e-2A^(?-?o) ' 



6 Comparison with Frohlich's Method II: The 
Degenerate Case 

If systems with many electrons and phonons are considered in general it will be 
not possible to avoid that the energy difference a^^q becomes zero. More precisely 
there will exist electron and phonon momenta satisfying the relation 

kq cos 6 



m 



± Vsq = 0, 



where 6 denotes the angle between k and q and Vs the velocity of sound. The 
Frohlich-transformation is not defined for these k, q. Thus the corresponding 
couplings Mk,q are not eliminated in this approach. 

But even in these cases the elimination described by the flow equations works. 
To verify this the asymptotic behavior of ujq has to be taken into account. 
For ak,q the following asymptotic expansion holds: 

afc,g(0 ~ «fc,<7(00) + Cg ■ / > 1. 



Thus for ttfc g(oo) = one obtains 



M.,^C-[^) =C.^-j , (49) 

with an appropriate constant C. In this case the decay is rather algebraic than 
exponential. 

Hence it is not important that ak^q becomes zero for some special I. It is only 
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important that it is not zero for all /. Thus the elimination of M^^q is a direct 
consequence of the renormalization of uOq. 

A similar argumentation was used by Kehrein, Mielke and Neu |P for the spin- 
boson model. There the authors argued that the coupling to the bosonic bath 
always is eliminated because of the renormalization of the tunneling frequency. 

A final remark shall be made on the remaining singularities of the electron- 
electron interaction ([T9|) . By taking the renormalization of ujq into account only 
the special values ^^^^(oo) = /9fe'^_q(oo) = are critical. It is easily seen that in 
this case 

Vk,k',q{l) = Vk,k',q{k) + D -logl/lQ, D := const. (50) 

holds for / ^ 1. Thus there will remain some singularities in our approach. But 
as long as these special values are multiplied with functions which stay finite for 
/ oo they will have no influence^. 

7 Conclusions 

It has been shown that it is possible to eliminate the electron-phonon interac- 
tion by using continuous unitary transformations. This elimination causes the 
renormalization of the coupling functions of the Hamiltonian described by the 
flow equations. This change of the couplings corresponds to the renormalization 
of the one-particle energies and to the generation of an effective electron-electron 
interaction. In order to set up these differential equation the generated new 
interactions have to be neglected. 

By analyzing these flow equations some different approaches can be chosen. 
First it is possible to apply the approximation used by Frohlich. In this case 
terms of higher order than |M(0)p can be neglected. Furthermore the possibil- 
ity of degeneracies has to be excluded. By doing so the flow equations can be 
solved exactly. By comparing the results with Frohlich's one sees that the induced 
electron-electron interactions differ. The origin of this difference lies in the fact 
^This will typically happen if one calculates averages of observables. 
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that in our approach all couplings depend on /. Thus by going over to a formula- 
tion similar to Frohlich's, i.e. carrying out the /-ordering of t], it can be seen that 
the flow of the coefficients changes the generator of the unitary transformation. 
Our result looks much friendlier than Frohlich's interaction. But nevertheless it 
is not clear by now in which sense our result is 'better' or even 'the best' one can 
obtain. 

It is worth mentioning that Kehrein and Mielke obtain similar modifications 
by eliminating the hybridization term in the single impurity Anderson model by 
continuous unitary transformations [0]. The authors show that their approach 
generates a spin-spin interaction which differs from the one obtained by the fa- 
mous Schrieffer- Wolff transformation. Their induced interaction is also less sin- 
gular. But even more important their result shows the right high-energy cutoff. 
Thus at least for this model a physical criterion exists to decide which interaction 
is 'better'. 

Besides this direct comparison with Frohlich's results it is possible to prove 
that our approach indeed eliminates the interactions between one phonon and 
one electron. For this purpose the influence of the renormalization of Ug on the 
elimination of Mk,q has to be taken into account. By doing so the renormalization 
of ujq in the asymptotic regime is described by an integro-differential equation. 
This equation can be transformed into a differential-equation which is much easier 
to analyze. By using its flrst integral the behavior of Ug for large / can be 
characterized. As a consequence of the properties of Ug the electron-phonon 
coupling is always eliminated even in the case of degeneracies. Furthermore it is 
shown that the remaining singularities in the electron-electron interaction are of 
a very weak nature. 

Acknowledgment. The authors would like to thank S. Kehrein and A. 
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